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Here we report the measurement of water vapor in Titan's stratosphere using the Cassini Composite 
Infrared Spectrometer (CIRS, Flasar, F.M. et al. [2004], Space Sci. Rev. 115, 169-297). CIRS senses water 
emissions in the far infrared spectral region near 50 pm, which we have modeled using two independent 
radiative transfer codes (NEMESIS (Irwin, P.G.J. et al. [2008], J. Quant. Spectrosc. Radiat. Trans. 109, 1136- 
1150) and ART (Coustenis, A. et al. [2007], Icarus 189, 35-62; Coustenis, A. et al. [2010], Icarus 207, 461- 
476). From the analysis of nadir spectra we have derived a mixing ratio of 0.14 ± 0.05 ppb at an altitude of 
97 km, which corresponds to an integrated (from 0 to 600 km) surface normalized column abundance of 
3.7 ± 1.3 x 10 14 molecules/cm 2 . In the latitude range 80°S to 30°N we see no evidence for latitudinal vari- 
ations in these abundances within the error bars. Using limb observations, we obtained mixing ratios of 
0.1 3 ± 0.04 ppb at an altitude of 1 1 5 km and 0.45 ± 0.1 5 ppb at an altitude of 230 km, confirming that the 
water abundance has a positive vertical gradient as predicted by photochemical models (e.g. Lara, L.M., 
Lellouch, F., Lopez-Moreno, J.J., Rodrigo, R. [1996], J. Geophys. Res. 101(23), 261; Wilson, E.H., Atreya, 
S.K. [2004], J. Geophys. Res. 109, E6; Horst, S.M., Vuitton, V., Yelle, R.V. [2008], J. Geophys. Res., 113, 
E10). We have also fitted our data using scaling factors of ~0.1-0.6 to these photochemical model pro- 
files, indicating that the models over-predict the water abundance in Titan’s lower stratosphere. 

© 2012 Elsevier Inc. All rights reserved. 


1. Introduction 

Water is present in its various forms in many regions of the So- 
lar System, from the atmospheres of the inner planets and shadows 
of lunar craters, to the mantles of icy satellites and beyond to the 
Kuiper Belt Objects (KBOs) and Oort Cloud Comets. Liquid water 
is also an essential ingredient for life on Earth and a potential clue 
in the search for life or habitability conditions in the rocks of Mars, 
the internal ocean of Europa or Titan, and the volcanic vents of Enc- 
eladus. On Titan, Saturn’s largest satellite that hosts a dense nitro- 
gen-dominated atmosphere, water is a trace species in the 
atmosphere. However, water plays a significant role since it is 
one of the sources of oxygen for the observed active photochemis- 
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try on Titan (e.g. Lara et al., 1996; Wilson and Atreya, 2004; Horst 
et al., 2008). 

Titan’s known oxygen compounds to date are carbon monoxide 
(CO, ~47 ppm), carbon dioxide (C0 2 , ~15 ppb) and water vapor 
(H 2 0), where the abundances are quoted for the low-latitude 
stratosphere (de Kok et al., 2007a). C0 2 was first detected by Voy- 
ager 1 (Samuelson et al., 1983), while CO was first seen by ground- 
based observations in the near-lR (Lutz et al., 1983). Subsequent 
observations in the sub-millimeter led to controversy as to 
whether CO was well-mixed or not (Hidayat et al., 1998; Gurwell, 
2004). CO emission lines were later observed by the Cassini Com- 
posite Infrared Spectrometer (CIRS), thus improving the previous 
abundance estimate (de Kok et al., 2007a; Teanby et al„ 2009). 
Water was first detected in Titan’s atmosphere by the Infrared 
Space Observatory (ISO) in 1997 (Coustenis et al., 1998). Two lines 
near 40-pm observed by the Short Wavelength Spectrometer 
(SWS) were modeled using a uniform mixing ratio above the con- 
densation level and a value of 0.4 ppb was detected. An early 
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attempt to measure H 2 0 with Cassini CIRS was unsuccessful due to 
poor signal-to-noise (S/N) ratios in early versions of the calibration 
pipeline spectra and a limited number of available spectra. There- 
fore, only an upper limit of 0.9 ppb could be retrieved (de Kok et al., 
2007a). Since then, water emission in the CIRS data has been defi- 
nitely observed, albeit without deriving any further information on 
its abundance and distribution (Bjoraker et al., 2008). The Cassini 
Ion and Neutral Mass Spectrometer (INMS) detected H 2 0 in the 
upper atmosphere - between 950 and 1200 km - with a mixing ra- 
tio in the range of ~(0.4-3.4) x 10~ 5 (Cui et al., 2009). 

While the presence of these oxygen compounds is now well- 
established, some details about their origin remain to be deter- 
mined. Early photochemical models assumed that CO originated 
from episodic outgassing from Titan’s interior along with nitrogen 
(N 2 ) or ammonia (NH 3 ) and methane (CH 4 ), whereas water mole- 
cules entered the top of the atmosphere and were photochemically 
dissociated to produce hydroxyl radicals (OH) (Wong et al., 2002; 
Wilson and Atreya, 2004). The combination of OH and CO led to 
the production of C0 2 . However, Horst et al. (2008) have recently 
challenged this model, arguing instead that both CO and C0 2 are 
the result of upper-atmospheric chemistry that occurs between 
in-falling oxygen species reacting with carbon produced by CH 4 
photodissociation. In this hypothesis, water enters Titan's atmo- 
sphere either in the form of H 2 0 or OH (since the latter is quickly 
converted to H 2 0 within the atmosphere) together with oxygen (O 
and 0 + ). These forms of oxygen are thought to be deposited on Ti- 
tan at two different altitudes. 0 + ions have been observed flowing 
into Titan’s atmosphere (Hartle et al., 2006a, b) and they are 
thought to be deposited in the upper atmosphere around 
1100 km (Horst et al„ 2008) where their interaction with methyl 
(CH 3 ) radicals leads to the formation of CO. Water is instead depos- 
ited at 750 km due to micrometeoritic ablation (English et al., 
1996) where it is photolyzed to OH. The latter finally combines 
with CO to form C0 2 and possibly other complex species. 

Saturn’s rings and the icy satellites that surround the giant plan- 
ets, and also interplanetary dust, are probable sources of the water 
(oxygen) in Titan’s atmosphere and recent results from INMS 
indicate that the plumes of Enceladus are the dominant source 
(e.g. Dougherty et al., 2006). Observations and models of the neutral 
H 2 0, OH, and O torus formed from the Enceladus plume show that 
material from Enceladus extends well beyond Titan’s orbit (Melin 
et al., 2009; Cassidy and Johnson, 2010; Fleshman et al., 2012). Based 
on Herschel measurements of the Enceladus torus combined with 
modeling of the fate of the species within the torus, Hartogh et al. 
(2011) showed that the flux of 0/0 + into Titan is consistent with 
an Enceladus source for the oxygen seen in Titan CO, except for the 
fact that Enceladus does not seem to provide enough OH/H 2 0. 

In this paper, we analyze the spectra acquired by CIRS in the far 
infrared spectral region in order to retrieve the water vapor vertical 
or spatial distribution in Titan’s atmosphere. CIRS has been acquir- 
ing spectra of Titan since the beginning of the Cassini prime mission 
(July 2004). After 2 years of the extended mission (XM), which in- 
cluded the 2009 equinox, in July 2010 Cassini entered in the Solstice 
Mission (SM), which is scheduled to last until 2017. Since an upper 
limit for H 2 0 was reported by de Kok et al. (2007a) there has been a 
considerable increase of the number of data collected by CIRS and 
significant improvements to their calibration. The increased signal 
to noise (S/N) ratio not only permits a definitive detection of H 2 0 
from the analysis of CIRS far infrared spectra, but it allows us to con- 
strain its vertical and latitudinal profile. 

2. Selected dataset 

CIRS (Flasar et al., 2004) is comprised of three Focal Planes 
observing in the spectral range 10-1400 cm 1 with spectral resolu- 


tions from 0.5 to 15.5 cm -1 . Focal Plane 1 detector (FP1) is charac- 
terized by a circular field of view (FOW) of 3.9 mrad. It records data 
in the far infrared spectral range (10-600 cm -1 ) with a spectral 
resolution of 0.5 cm -1 , allowing us to observe the water vapor sig- 
nature, and by modeling, to retrieve its abundance. Water presents 
its rotational lines in the CIRS FP1 spectral region up to 400 citT 1 , 
with the strongest and most visible lines in the range positioned 
between 90 and 260 cm -1 . We focus here on the range from 150 
to 260 cirr 1 for the water detection, as this is the range of maxi- 
mum responsivity of FP1. At lower wavenumbers the on-board 
electronics of CIRS create a moving interference spike that can af- 
fect the spectrum up to 150 cm -1 . Therefore, we exclude wave- 
numbers shorter of 150 cm -1 . We use data from two different 
types of observations to obtain independent measurements: the 
far infrared on-disk integrations (FIRNADCMP) and the far infrared 
limb integrations (FIRLMBINT). Water is a trace species with rela- 
tively weak lines and therefore it cannot be observed in an individ- 
ual spectrum. An average of a few thousand spectra of on-disk 
observations and a few hundred spectra of limb observations is 
necessary to achieve sufficient signal-to-noise (Fig. 1). 

Limb observations have the FP1 Focal Plane centered around 
two different altitudes - hereafter limb 1 and 2 - and are therefore 
used to constrain the water vapor abundance in the stratosphere 
around 115 and 230 km respectively, well above the 45 km tropo- 
pause. Since the contribution functions of water for on-disk obser- 
vations peak around 97 km (Fig. 2), the retrieved water vapor 
abundance derived from these measurements can be compared 
with the lowest altitude targeted by our limb integrations around 
115 km. 

For the water detection and retrieval of quantitative informa- 
tion together with possible latitudinal variations, multiple Titan 
flybys must be utilized to enhance the signal. To date, 35 limb inte- 
grations of approximately 1 h in duration (~60 high-resolution 
spectra) have been obtained covering latitudes from 87°S to 
80°N. The nadir integrations are more numerous (about 92 suc- 
cessfully executed, of typical duration 5 h, ~300 spectra) as they 
occur in a less contested observing time further from the Titan 
closest-approach period. They also have a more or less complete 
spatial coverage of Titan’s latitudes and longitudes with an average 
footprint size of ~15° great circle arc. 

We focus on one season of on-disk observations acquired from 
December 2004 to December 2008 (northern winter on Titan) in 
order to reach a compromise between obtaining a large number 
of spectra and a sufficiently homogeneous dataset. Inside this time 
period for on-disk observations acquired from a maximum dis- 
tance of 300,000 km and with a maximum emission angle of 60°, 
we selected latitudinal bins (80-45°S, 45-10°S, and 0-30°N) cen- 
tered around three latitudes for which observation-derived tem- 
perature profiles were available (see model description in Section 
3 and Fig. 3). The number of spectra averaged in these latitudinal 
bins were approximately 1700, 3800 and 7000 respectively and 
their average emission angles were 35°, 38° and 34°. 

CIRS limb spectra are acquired in much smaller numbers, there- 
fore to reach a sufficient signal-to-noise ratio we consider only one 
average of about 320 spectra acquired from December 2004 and 
September 2009, encompassing the entire southern hemisphere 
and mid-latitudes within the range 90°S to 20°N. During this time 
period and at these latitudes data can be considered quite homoge- 
neous as shown in Teanby et al. (2010), therefore we model this 
average using a temperature profile retrieved for 1 5°S. We exclude 
from this analysis the higher northern latitudes where the strato- 
spheric temperature profile changes significantly. We have also se- 
lected data acquired by the spacecraft at a range less than 
45,000 km in order to limit the size of the projected detector foot- 
print on the limb to less than 1 50 km. 
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Fig. 1. (a) In black, the average of CIRS far-IR on-disk observations is plotted (~7000 spectra acquired from December 2004 to December 2008 in the latitudinal range of 0- 
0°N), limb observations centered around 115 and 230 km (respectively ~320 and ~280 spectra acquired from December 2004 to September 2009 in the latitudinal range of 
90-20°S) and their fit (in green, blue, red respectively) assuming a constant water mole fraction above the condensation altitude. In (b) the retrieval spectral range is shown 
with the main water lines indicated by vertical dotted lines. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this 
article.) 


Normalised Water Contribution Function 
0.0 0.2 0.4 0.6 0.8 1.0 



Fig. 2. Contribution functions of the different atmospheric layers to the water vapor line emission computed for four wave numbers. The solid line shows a temperature 
profile for 15°N. 


3. Data analysis and model 

In the selected part of the FP1 spectral range, Titan’s spectrum is 
formed by (i) the contribution of thermal emission of the surface 


and atmospheric layers, (ii) the seven pairs (Anderson and Samuel- 
son, 2011) of collision induced absorption (CIA) opacities between 
the main atmospheric molecules - N 2 , CH 4 and H 2 - due to Titan’s 
dense lower atmosphere, (iii) the photochemical aerosol plus 
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Fig. 3. Atmospheric vertical temperature-pressure profiles retrieved from CIRS data for three latitudes (58°S, 15°S, 15°N) from Anderson and Samuelson (2011) from the 
surface to 3.3 x 1 0 7 bar, corresponding to an altitude range of 0-600 km. A priori constant with altitude water vapor profiles are also shown for the three temperatures 
profiles. The constant mixing ratio profile along the atmosphere (a priori assumed to be 0.1 ppb) is decreased to follow the saturation vapor pressure curve in the lower 
stratosphere below the condensation altitude ('cond. alt.’ in the figure), where the relative humidity reaches 100%. Using the saturation vapor pressure equation of water over 
ice of Murphy and Koop (2005) and assuming the temperature profiles retrieved at 58°S, 15°S, 15°N we find a condensation altitude of about 104, 97 and 93 km respectively. 


stratospheric condensates, and (iv) the ro-vibrational emission 
lines of atmospheric species seen by CIRS at the latitudes consid- 
ered in our study: CH 4 , CO, H 2 0, C 4 H 2 . 

These quantities were used as input to the NEMESIS retrieval 
code (Irwin et al., 2008) to perform a combination of forward mod- 
el computation and retrieval scheme based on the method of opti- 
mal estimation (Rodgers, 2000). The computation of the forward 
model spectrum used the correlated-k approach of Lacis and Oinas 
(1991) and included a Hamming apodization of Full Width at Half 
Maximum (FWHM) of 0.5 cm -1 to reproduce the instrumental line 
shape. The retrieval scheme was used to optimize the fits and 
determine the model free parameters including the water vapor 
abundance. This method was successfully applied to model the 
FP1 spectrum in Cottini et al. (2012) to retrieve surface tempera- 
ture. We solve the radiative transfer equation for 147 spherical 
atmospheric layers, using as a source function the thermal emis- 
sion of the surface, for which a unit surface emissivity is assumed, 
and that of the atmospheric layers. The retrieval algorithm then 
iteratively computes a synthetic spectrum, compares it to the data 
and after applying a cost function, determines the best estimate for 
the physical parameters in the model - the stratospheric aerosol 
profile and any necessary adjustments to the temperature profile 
and the mole fraction of included atmospheric gases. The cost func- 
tion includes two components: one that measures the quality of 
the fit to the spectra (similar to a % 2 test) and another that mea- 
sures the deviation of the retrieved parameters from a set of a pri- 
ori quantities. 

The continuum due to the CIA was calculated according to Bory- 
sow and Frommhold (1986a,b,c, 1987), Borysow (1991), and Bory- 
sow and Tang (1993). For the N 2 -CH 4 pair, we used CIA coefficient 
values increased by 50% as required to fit the continuum of the 
Cassini Descent Imager Spectral Radiometer (DISR) data (Tomasko 
et al., 2008) and the CIRS spectra (de Kok et al., 2010). 

We have modeled the haze emission/absorption using the 
extinction cross sections of the hazes included in de Kok et al. 
(2007b). Since scattering is negligible at these wavelengths for par- 
ticles smaller than few microns, we have omitted it from our 
computations. 

We have adopted the atmospheric vertical temperature-pres- 
sure profiles retrieved from CIRS data for three latitudes (58°S, 
15°S, 15°N) from Anderson and Samuelson (2011) from the surface 
to 3.3 x 10 -7 bar, corresponding to an altitude range of 0-600 km 


(Fig. 3). Spectroscopic information for the gas rotational lines in the 
far-infrared range was extracted from the HITRAN 2004 database 
(Rothman et al., 2005). For CH 4 we have adopted the revised mole 
fraction of 1.48% in the stratosphere (Niemann et al., 2010) ac- 
quired by the Gas Chromatograph Mass Spectrometer (GCMS) on 
the Huygens probe during in its descent to Titan’s surface. In the 
stratosphere for H 2 we assumed a uniform volume mixing ratio 
of 0.1% (Courtin et al., 2008). The geometry of the observations 
was also included in the computations. 

An accurate model of the instrumental FOV is required to suc- 
cessfully reproduce CIRS spectra and measure water abundance. 
The FP1 FOV is circular and has a sensitivity with a quasi-exponen- 
tial decrease from the center to the edge and a FWHM of 2.4 mrad. 
For on-disk spectra the homogeneity of the field of view usually 
permits simple modeling with a single ray calculated for the detec- 
tor center. For limb spectra, we have to take into account the rapid 
decrease in atmospheric density with height and the variations of 
temperature and gas volume mixing ratio profile with altitude. In 
such cases, the FOV is not assumed to be uniform and a multiple 
ray model is required to fit the data. We modeled the FOV using 
the minimum number of rays for which the synthetic spectrum 
and the water retrieval computation results became stable; this 
corresponds to nine rays with a step in altitude of 25 km. We also 
recomputed some of the results using 39 rays (step of 5 km) in or- 
der to show a smoother limb contribution function. The spectral 
radiance measured by the FP1 detector is modeled by a convolu- 
tion of the emerging radiance at each point in the FOV (as de- 
scribed in Nixon et al. (2009) and in Teanby and Irwin (2007)), 
weighted by a response function for CIRS FP1 detector. This beam 
profile was determined for CIRS FP1 (Flasar et al., 2004) using Jupi- 
ter as a point source. As 95% of the integrated response is contained 
in a radius of 1 .95 mrad from the FOV center, the detector observes 
a maximum altitude range of about 70 km. 

Independent line-by-line calculations to simulate the same on- 
disk FP1 selections were also made using the Atmospheric Radia- 
tive Transfer (ART) code that has recently been applied to CIRS data 
in Coustenis et al. (2010, 2007). The code uses the most recent aer- 
osol extinction dependence inferred from Vinatier et al. (2012) and 
temperature profiles derived by fitting the v 4 methane band at 
1304 cm -1 in FP4 averages taken at similar conditions as the FP1 
spectra. The spectroscopic parameters for all the observed mole- 
cules and isotopes are from GE1SA 2009 (Jacquinet-Husson et al., 
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2011) and HITRAN 2008 (Rothman et al., 2009). The results from 
these two different codes are very similar and their difference is 
smaller than the error bars on the data. 

In Fig. la the on-disk and two limb observation averages are 
shown together with their fits. Fig. lb shows only the spectral 
range used for the water line analysis. 

4. Results 

We have retrieved the water vapor abundance from both on- 
disk and limb data assuming different vertical profiles. These were 
a constant water mixing ratio profile and three vertically increas- 
ing profiles from recent photochemical models: (a) Horst et al. 
(2008), (b) Wilson and Atreya (2004) and (c) Lara et al. (1996). 
The latter profile was adopted in Coustenis et al. (1998) for the first 
water detection on Titan by ISO. The constant mixing ratio profile 
(a priori assumed to be 0.1 ppb) was decreased to follow the satu- 
ration vapor pressure curve in the lower stratosphere below the 
altitude where the relative humidity reaches 100% (see Fig. 3). 
Using the saturation vapor pressure equation of water over ice of 
Murphy and Koop (2005) and assuming the temperature profile re- 
trieved at 15°N we find a condensation altitude of 93 km. For tem- 
perature profiles retrieved at 15°S and 58°S we find respectively 
condensation altitudes of 97 km and 104 km. We have also com- 
puted the contribution functions - normalized inversion kernels 
- showing the sensitivity of each atmospheric layer to a variation 
of the H 2 0 mixing ratio. These contribution functions were com- 
puted for each profile and for all of the most intense water lines 
in order to provide an altitude range of validity of the retrieved val- 
ues (Figs. 2 and 4). Fig. 2 shows the contribution functions for on- 
disk observations computed at four different wave numbers; at 
254 cm 1 (one of the two lines used for the ISO water retrieval) 
the upper shoulder of the contribution function is wider and sensi- 
tive to higher altitudes compared to the other wave numbers used 
for the water retrieval in this work. In our case the fit of the 
254 cm 1 line improves when using a profile increasing with alti- 
tude rather than a constant profile. In Fig. 4 we show only the con- 
tribution functions computed at wavenumber 202.75 cnr 1 , where 


the most intense water line in the CIRS spectrum occurs (discount- 
ing the line at 1 50.5 cnr 1 that is unusable due to an instrumental 
interference). 

For water retrievals obtained using a scaled constant water 
profile we show the retrieved mixing ratio values at the altitude 
where the water contribution function peaks for the assumed 
profile (Table 1 ). For the vertical error we use the Full Width Half 
Maximum (FWHM) of the contribution function for the corre- 
sponding water profile. We also retrieve a scaling factor to the 
water profile from each of the photochemical models considered 
in this work (Table 1 ). 

4.1. On-disk water retrieval 

To measure the water abundance from the on-disk average (0- 
30°N) data we first use a constant water profile. We retrieve a vol- 
ume mixing ratio of 0.14 ± 0.05 ppb at 97 km (FWHM 93-130 km). 
This value corresponds to a surface-normalized H 2 0 total column 
density of 3.7 ± 1.3 x 10 14 molecules/cm 2 . 

The largest source of error for the on-disk observations is due to 
the fact that the maximum of the contribution function occurs in 
the region where the water abundance is rapidly changing due to 
condensation. Other sources of error include small variations of 
the temperature profile in the stratosphere, random noise from 
the detectors, and a small dependence in altitude sensitivity with 
wavenumber. 

We also fit the water lines for the three other water vertical dis- 
tribution profiles (Lara et al., 1996; Wilson and Atreya, 2004; Horst 
et al., 2008) and obtained the necessary scale factors to fit the data, 
which are shown in Table 1. These values, ranging between 0.11 
and 0.63 times the considered profiles, show the retrieved water 
mole fraction to be less than predicted from these previous models. 

We have analyzed two additional latitudinal bins from 45°S to 
10°S and from 80°S to 45°S, centered on the latitudes correspond- 
ing to the temperature profiles previously retrieved from CIRS at 
15°S and 58°S respectively (Anderson and Samuelson, 2011). See 
Fig. 3. The observed water mixing ratio indicates the absence of 
any significant latitudinal variations within the error bars in these 



Normalised Water Contribution Function 



Normalised Water Contribution Function 


Normalised Water Contribution Function 




Fig. 4. Contribution functions of water vapor line emission and temperature profile for 15°N. In (a-d) we show the contribution functions computed at 202.75 cm -1 for 
different water profile models: (a) a constant water vertical profile; (b) profile from model D in Horst et al. (2008); (c) profile from Wilson and Atreya (2004) and (d) profile 
from Lara et al. (1996). 
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Table 1 

Water vapor abundance results. 


Water vapor retrieved mole fractions 

On-disk average (0-30°)N 

Limb retrieval 1 

Limb retrieval 2 

Constant VMR profile 

Scaling factor to H 2 0 profile from Horst (model D) 
Scaling factor to H 2 0 profile from Wilson-Atreya 
Scaling factor to H 2 0 profile from Lara 

(0.14 ± 0.05) ppb at 97±f km 
0.18 ±0.05 at 1 1 8 km 
0.14 ±0.05 at 118'“ km 
0.48 ±0.07 at 115f^ km 

(0.13 ±0.04) ppb at 1 1 5+1° km 
0.14 ±0.05 at 129 ?) km 
0.13 ±0.05 at 129f?4 km 
0.63 ±0.07 at 133f?) km 

(0.45 ± 0.15) ppb at 230^ km 
0.23 ± 0.07 at 232fg km 
0.18 ±0.08 at 222^3 km 
0.45 ± 0.08 at 247+^ 0 km 



Fig. 5. Water vapor mole fraction retrieved from CIRS on-disk and two limb 
observations assuming a water profile constant with altitude over the condensation 
level. Water profiles from three photochemistry models are also shown for 
comparison. The Horst et al. (2008) water vapor profile was derived assuming 
two different eddy diffusion coefficients - 100 cm 2 s 1 and 400 cm 2 s 1 (the second 
being the one recommended in their model; dotted curve plus solid curve). Also the 
profiles in Wilson and Atreya (2004; dashed curve) and in Lara et al. (1996; dot- 
dash curve). 

latitude ranges. It should be stressed that in this work we did not 
analyze the water stratospheric content at high northern latitudes 
that were experiencing winter during this time period. To model 
spectra at these latitudes for a large average is particularly com- 
plex since the stratospheric temperature was changing quickly 
with latitude. In addition the temperature profiles for high north- 
ern latitudes are not yet available for the lower stratospheric re- 
gion sensed by the CIRS water lines. 

The results obtained simultaneously for on-disk data using the 
independent line-by-line ART code simulations of the same FP1 
selections confirm the retrieved water vapor value reported above 
within the error bars. 

4.2. Limb water retrieval 

The measurement of water vapor obtained by modeling the 
limb 1 spectrum under the assumption of a constant mixing ratio 
profile is equal to 0.13 ± 0.04 ppb. According to the position of 
the peak of the corresponding contribution function, the radiance 
mostly originates from a region centered at 115 km (FWHM 95- 
165 km). Modeling the limb 2 spectrum we retrieved a water mix- 
ing ratio using a constant water profile of 0.45 ± 0.1 5 ppb at an alti- 
tude of 230 km (FWHM 190-275 km). These values indicate an 
increase of the water mole fraction with altitude in the strato- 
sphere from 115 km to 230 km of about three times. 

The scaling factors to the model water profiles obtained from 
the two types of limb retrievals are shown in Table 1 and again 
illustrate the smaller amount of stratospheric water vapor detected 
by CIRS with respect to the ones predicted by the models consid- 
ered in this study. 

5. Conclusions 

In this work we modeled CIRS data using a constant-with- 
height water vapor profile and assigned the retrieved mixing ratio 
to the altitude where the contribution function peaks. 


By combining on-disk and limb observations we are able to con- 
strain the vertical profile of water in the stratosphere from 12 mbar 
to 1 CT 3 mbar, corresponding to altitudes between 93 and 280 km 
(considering the widths of the contribution functions). 

In Fig. 5 we compare our water vapor retrieved values with the 
models of (Lara et al., 1996; Wilson and Atreya, 2004; Horst et al., 
2008). 

The measurement of the stratospheric vertical profile of water 
adds useful constraints to the photochemical models of Titan’s 
atmosphere. Qualitatively, the increase of the water mixing ratio 
with altitude is in agreement with an external source of oxygen 
and a lower altitude sink due to condensation. However, quantita- 
tively, our retrieved abundance seems to be less (from ~0.1 to 
~0.6) than predicted from the models considered in this work 
(see Table 1). We also observe that since the scaling factors to 
the photochemical models (a)-(c) in Table 1 are slightly different 
for the two limb altitudes, it implies that these models might have 
a slope for H 2 0 not quite consistent with CIRS data. However, due 
to the rapid variation of the water vapor mole fraction with alti- 
tude in the atmospheric region where water freezes, and where 
CIRS is actually observing (on-disk and limb 1 spectra), we should 
be cautious in assigning a slope to the water profile. 

Coustenis et al. (1998) fitted the ISO data with a scaling factor to 
the Lara et al. (1996) water profile of 0.4^2- We fitted the CIRS 
data assuming the same profile multiplied by a scaling factor of 
0.48 ± 0.07 (corresponding to a water column density of 
3.8 ± 1.0 x 10 14 molecules cm -2 ). This result agrees with the scal- 
ing factor from the ISO analysis although this agreement may be 
fortuitous because most of the ISO emission originates from above 
300 km, with the contribution functions peaking around 400 km 
(Fig. 3 of Coustenis et al. (1998)). Since the beam size of ISO was 
much larger than Titan there is a strong emission from the limb 
occurring at high altitudes where 0.4 times the Lara et al. (1996) 
profile is used. Due to its higher spatial resolution, the contribution 
functions for the CIRS nadir selection using the Lara et al. (1996) 
profile cover the range 95-145 km (at half maximum). 

ISO also retrieved a water vapor abundance of 0.4 ppb assuming 
a constant mole fraction above the condensation level (Coustenis 
et al., 1998); we assign the same relative error bars as those de- 
rived by ISO for the scaled Lara et al. (1996) profile: +0.3 and 
-0.2. From CIRS on-disk observations we retrieved a volume mix- 
ing ratio of 0.14 ± 0.05 ppb around 97 km for latitudes 0-30°N. Our 
retrieval is only marginally consistent with the ISO determination 
of 0.4 4 q 2 ppb above the water vapor condensation altitude. How- 
ever it is necessary to make the assumption that the ISO results 
pertain to the same condensation altitude that we retrieve because 
unfortunately it is not reported in the ISO paper. 

We now discuss how our retrieved water abundances compare 
with current models of oxygen photochemistry on Titan. The re- 
trieved scaling factors are all less than one, implying that there 
are still uncertainties in our understanding of oxygen processes 
on Titan. 

In photochemical models of Titan prior to 2000 (e.g. Lara et al„ 
1996), in order to allow atmospheric production of CO, it was pos- 
tulated that CO could be produced through a chemical reaction be- 
tween OH (available from H 2 0 influx into the upper atmosphere) 
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and CH 3 . However, it was emphasized by Wong et al. (2002) that 
this reaction does not produce CO as previously assumed but in- 
stead it produces H 2 0 (Pereira et al., 1997). 

Hence an influx of H 2 0 or OH does not produce any significant 
abundance of CO and therefore C0 2 can be produced by an H 2 0 in- 
flux only with CO already present (OH + CO -> C0 2 + H). For this 
reason these models were unable to reproduce the observed CO 
abundance and were substituted by other models that suggest 
the existence of primordial CO in the atmosphere (Wilson and 
Atreya, 2004), or instead CO produced in the atmosphere using 
an external influx of 0 + rather than H 2 0 or OH (Horst et al., 2008). 

In the pre-Cassini model of Wilson and Atreya (2004) water is 
photolyzed to OH, which combines with CO to form C0 2 and other 
complex species. In this model, CO is assumed to be primordial on 
Titan and the water abundance profile derives from the amount nec- 
essary to form the observed C0 2 . This assumption was challenged by 
the Horst et al. (2008) model, in which oxygen species are assumed 
to arrive from outside the Moon and form carbon monoxide as well 
as carbon dioxide in the atmosphere. The values of the input fluxes 
of O and OH were adjusted to reproduce the observed abundances of 
CO and C0 2 . In the Horst et al. (2008) model water profiles were 
produced for six different values (from K 0 = 100 cm 2 s -1 to 
/Co = 1000 cm 2 s _1 ) of the eddy coefficient in the lower atmosphere, 
since the stratospheric abundances of photochemically produced 
species are highly dependent on this parameter. As shown in 
Fig. 5, the water abundance retrieved in our study is best fit by the 
water profile with the lowest eddy diffusion coefficient value con- 
sidered in their model (K 0 =100 cm 2 s~'). This is lower than the va- 
lue (/( 0 = 400 cm 2 s -1 ) they identified as best reproducing C1RS 
observations of hydrocarbon species and adopted in this work for 
comparison with our retrievals. However our results show that even 
the Horst model with K 0 =100 cm 2 s _1 still has excessive water at 
the altitudes of our measurements. 

Therefore, our work clearly points towards further refinement of 
oxygen chemistry in photochemical models of Titan’s atmosphere. 
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